Beaver dam analogues drive heterogeneous groundwater–surface water interactions

Beaver dam analogues (BDAs) are a cost-effective stream restoration approach that leverages the recognized environmental benefits of natural beaver dams on channel stability and local hydrology. Although natural beaver dams are known to exert considerable influence on the hydrologic conditions of a stream system by mediating geomorphic processes, nutrient cycling, and groundwater–surface water interactions, the impacts of beaver-derived restoration methods on groundwater–surface water exchange are poorly characterized. To address this deficit, we monitored hyporheic exchange fluxes and streambed porewater biogeochemistry across a sequence of BDAs installed along a central Wyoming stream during the summer of 2019. Streambed fluxes were quantified by heat tracing methods and vertical hydraulic gradients. Biogeochemical activity was evaluated using major ion porewater chemistry and principal component analysis. Vertical fluxes of approximately 1.0 m/day were observed around the BDAs, as was the development of spatially heterogeneous zones of nitrate production, groundwater upwelling, and anaerobic reduction. Strong contrasts in hyporheic zone processes were observed across BDAs of differing sizes. This suggests that structures may function with size-dependent behaviour, only altering groundwater–surface water interactions after a threshold hydraulic step height is exceeded. Patterns of hyporheic exchange and biogeochemical cycling around the studied BDAs resemble those around natural beaver dams, suggesting that BDAs may provide comparable benefits to channel complexity and near-stream function over a 1-year period.     

Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

海洋表层沉积物中金属镉的来源解析

利用2015年8月份对长江口及其邻近海域表层沉积物的监测数据,基于主成分分析/绝对主成分分数(PCA/APCS)受体模型定量解析了重金属元素镉(Cd)的可能来源,并结合地统计学插值了Cd的源贡献量的空间分布状况,结果表明沉积物中镉污染主要存在3个可能来源,源头及贡献率分别为工业污染(18.8%)、陆地径流输入(66.0%)、生物活动等自然因素(13.6%),并且各个源头贡献量具有不同的空间分布状况,其中工业污染的高值区主要集中在靠近陆地区域,陆地径流输入的分布呈现由近岸向外海逐减降低的特征,生物活动等自然因素的高值区主要集中在远离陆地的外海区域。     

Ephemeral rollover points and clinothem evolution in the modern Po Delta based on repeated bathymetric surveys

Reconstructions of ancient delta systems rely typically on a two-dimensional (2D) view of prograding clinothems but may miss their three-dimensional (3D) stratigraphic complexity which can, instead, be best documented on modern delta systems by integrating high-resolution geophysical data, historical cartography, core data and geomorphological reconstructions offshore. We quantitatively compare three precisely positioned, high-resolution multi-beam bathymetry maps in the delta front and pro delta sectors (0.3 to 10 m water depth) of Po di Pila, the most active of the modern Po Delta five branches. By investigating the detailed morphology of the prograding modern Po Delta, we shed new light on the mechanisms that control the topset to foreset transition in clinothems and show the temporal and spatial complexity of a delta and its pro delta slope, under the impact of oceanographic processes. This study documents the ephemeral nature of the rollover point at the transition between sandy topset (fluvial, delta plain to mouth-bar) and muddy seaward-dipping foreset deposits advancing, in this case, in >20 m of water depth. Three multibeam surveys, acquired between 2013 and 2016, document the complexity in space and time of the topset and foreset regions and their related morphology, a diagnostic feature that could not be appreciated using solely 2D, even very high-resolution, seismic profiles. In addition, the comparison of bathymetric surveys gathered with one-year lapses shows the migration of subaqueous sand dunes on the clinothem topset, the formation of ephemeral cut-and-fill features at the rollover point (few m below mean sea level), the presence of collapse depressions derived by sagging of sediments and fluid expulsion (possibly induced by storm waves) on the foreset, and splays of sand likely reflecting gravity flows on the lower foreset. Though the modern Po Delta is anthropogenic in many respects, its subaqueous clinothem can be studied as a scale model for ancient clinothems that are less resolved geometrically and far less constrained chronologically.     

新疆和田地区耕地面源污染生态风险评价

随着化肥、农膜等在农业生产中的过量投入，耕地面源污染的程度随之加重。文章选取塔里木河流域上游和田地区为研究区域，依据P-S-R框架理论，构建和田地区耕地面源污染生态风险评价指标体系，加入土壤理化数据，使用生态风险评价模型对和田地区1980 年及2016 年耕地面源污染状况进行生态风险评价，运用耕地生态风险模型、生态风险转移矩阵、Arcgis分析和田地区耕地面源污染时空分异状况。研究结论如下：和田地区1980 年耕地生态风险等级均为II级或III级，呈“中间高，两侧低”分布；2016 年耕地生态风险等级上升至IV级或V级，呈“倒W型”分布，各县耕地面源污染程度较1980 年均有较大幅度的上升，其中墨玉县和于田县在2016 年耕地生态风险等级达到最高的V级，而民丰县因自身生态环境的强脆弱性，同样需要提高关注。根据面源污染“从源头治理”的原则，应切实推进和田地区耕地生态环境保护与治理，提高政府重视程度，增强技术指导，开展试点工作，改善和田地区耕地面源污染现状。     

徐闻南山镇红树林沉积物及红树植物中重金属的累积特征与迁移规律

为探究重金属在红树林沉积物及红树植物中的分布累积及迁移规律，选取了徐闻南山镇红树林为研究对象，通过测定红树林沉积物及红树植物不同部位（根、茎、叶）的重金属质量分数，运用富集因子、生物富集系数、转移系数及相关性分析等方法进行分析。结果表明：1）红树林沉积物重金属质量分数表现为铬（Cr）＞锌（Zn）＞镍（Ni）＞铜（Cu）＞铅（Pb）＞砷（As）＞汞（Hg）＞镉（Cd），为中等变异程度；除了镍（Ni）元素外，其余7种重金属未超过国家一级标准，除了铅（Pb）元素外，其余7种重金属均超过广东省土壤环境背景值，说明研究区沉积物中重金属具有一定的积累效应。2）沉积物中砷（As）、铜（Cu）、锌（Zn）、汞（Hg）、镍（Ni）、铬（Cr）富集因子值均＞1.5，说明受到轻微人为活动影响；各站位镍（Ni）富集因子值均＞5，结合研究区背景，反映了镍（Ni）受到自然和人为输入的共同影响。3）白骨壤体内重金属主要集中在根部，而红海榄体内重金属在根茎叶中分布相对均匀。白骨壤根茎叶部位的大多数重金属质量分数远高于红海榄，说明白骨壤对重金属的吸附能力比红海榄强。汞（Hg）集中分布在植物的叶片部位，且与其他重金属之间相关性不明显；推测汞（Hg）主要通过叶片吸收进入植物体内，与交通运输污染有关。4）不同红树植物对不同重金属富集能力各异，白骨壤对重金属的富集能力表现为：镉（Cd）＞砷（As）＞铜（Cu）＞锌（Zn）＞汞（Hg）＞铅（Pb）＞镍（Ni）＞铬（Cr），红海榄表现为：镉（Cd）＞铜（Cu）＞汞（Hg）＞锌（Zn）＞铅（Pb）＞砷（As）＞镍（Ni）＞铬（Cr）。白骨壤和红海榄对汞（Hg）的运移能力都较强；红海榄对镉（Cd）的富集能力和转运能力都较强，而白骨壤对镉（Cd）富集能力较强，转运能力却较弱，这说明红树植物对重金属元素的富集能力与转运能力不存在正比关系。     

基于ROMS和4DVAR的沿轨与网格化SSH数据同化效果评价

Remote sensing products are significant in the data assimilation of an ocean model. Considering the resolution and space coverage of different remote sensing data, two types of sea surface height(SSH) product are employed in the assimilation, including the gridded products from AVISO and the original along-track observations used in the generation. To explore their impact on the assimilation results, an experiment focus on the South China Sea(SCS) is conducted based on the Regional Ocean Modeling System(ROMS) and the four-dimensional variational data assimilation(4 DVAR) technology. The comparison with EN4 data set and Argo profile indicates that, the along-track SSH assimilation result presents to be more accurate than the gridded SSH assimilation, because some noises may have been introduced in the merging process. Moreover, the mesoscale eddy detection capability of the assimilation results is analyzed by a vector geometry–based algorithm. It is verified that, the assimilation of the gridded SSH shows superiority in describing the eddy's characteristics, since the complete structure of the ocean surface has been reconstructed by the original data merging.     

苏北浅滩沉积物中的大型绿藻微观繁殖体的垂直分布和物种多样性

本研究通过在2017年秋冬季组织实施的2个绿潮种源调查航次,研究苏北浅滩紫菜养殖区沉积物中的大型绿藻微观繁殖体的垂直分布和物种多样性。对沉积物中的微观繁殖体的定量实验结果表明:表层泥样中绿藻微观繁殖体的丰度最高,随着深度的增加迅速降低。在9月份表层沉积物中丰度为9±2ind./g,而在11月份表层丰度达到98±25ind./g;与9月份相比, 11月份浅滩沉积物中绿藻繁殖体数量呈现出显著升高的趋势。通过对51个绿藻样品的分子鉴定发现:沉积物中存在石莼属(Ulva)、尾孢藻属(Urospora)、盘苔属(Blidingia)等大型绿藻的微观繁殖体;数量占优是石莼属繁殖体,主要包括曲浒苔(Ulva flexuosa)、浒苔(Ulva prolifera)、扁浒苔(Ulva compressa)、缘管浒苔(Ulva linza);其中曲浒苔占据的比例最高,在9和11月份分别达到80.94%和73.33%。在2个调查航次中,都从沉积物中发现了绿潮浒苔的微观繁殖体,这些微观繁殖体是绿潮浒苔在苏北浅滩种源维持的重要方式和关键阶段,也构成了黄海浒苔绿潮连年暴发的"种子库"。     

Flexural deformation controls on Late Quaternary sediment dispersal in the Garo-Rajmahal Gap,NW Bengal Basin

Subsurface deformation is a driver for river path selection when deformation rates become comparable to the autogenic mobility rate of rivers. Here we combine geomorphology, soil and sediment facies analyses, and geophysical data of the Late Quaternary sediments of the central Garo-Rajmahal Gap in Northwest Bengal to link subsurface deformation with surface processes. We show variable sedimentation characteristics, from slow rates (<0.8 mm/year) in the Tista megafan at the foot of the Himalaya to nondeposition at the exposed surface of the Barind Tract to the south, enabling the development of mature soils. Combined subsidence in the Tista fan and uplift of the Barind Tract are consistent with a N-S flexural response of the Indian plate to loading of the Himalaya Mountains given a low value of elastic thickness (15–25 km). Provenance analysis based on bulk strontium concentration suggests a dispersal of sediment consistent with this flexural deformation—in particular the abandonment of the Barind Tract by a Pleistocene Brahmaputra River and the current extents of the Tista megafan lobes. Overall, these results highlight the control by deeply rooted deformation patterns on the routing of sediment by large rivers in foreland settings.     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.